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ABSTRACT 

Using long slit emission-line spectra we detect a fast outflow of ionized gas, with 
velocities up to 1000 km s _1 , in the nearby powerful radio galaxy 3C 293 (z — 0.045). 
The fast outflow is located about 1 kpc east of the nucleus, in a region of enhanced 
radio emission due to the presence of a distorted radio jet. We present results that 
indicate that this fast outflow is caused by a jet-ISM interaction. The kinematics of the 
outflowing ionized gas are very similar to those of a fast outflow of neutral hydrogen 
gas in this galaxy, suggesting that both outflows are the result of the same driving 
mechanism. While the mass of the outflowing ionized gas is about 1 x 10 5 M , the total 
Hi mass involved in the neutral outflow is about lOOx higher (10 7 M©). This shows 
that, despite the high energies that must be involved in driving the outflow, most of 
the gas remains, or becomes again, neutral. Other outflows of ionized gas, although 
not as pronounced as in the region of the enhanced radio emission, are also seen in 
various other regions along the axis of the inner radio jets. The regular kinematics of 
the emission-line gas along the major axis of the host galaxy reveal a rotating ionized 
gas disk 30 kpc in extent. 

Key words: galaxies: active - galaxies: individual: 3C 293 - galaxies: ISM - ISM: jets 
and outflows - ISM: kinematics and dynamics - line: profiles 



1 INTRODUCTION 

Nuclear activity can strongly influence the Interstellar 
Medium (ISM) in the centre of active galaxies. The strong 
radiation field from the central engine, as well as collimated 
jets of radio plasma, affect the ionization and kinematics 
of this ISM in a complex way. In addition, the situation 
can be further complica ted by the presence of st arburst- 
induced winds (e.g. lHeckman. Armus fc Milevll99d . and ref- 
erences therein), associated with young stellar populations 
that are known to be present in some radi o galaxies (e.g. 
lAretxaga et"al]|200li Iwills et alJl20ol .20041) . The feedback 
mechanisms of AGN and starburst activity are also impor- 
tant in the evolution of the host galaxies. AGN activity 
may, for example, regulate the correlation between the mass 
of t he central black h ole and the galaxy's bulge properties 
(e.g. lSilk fc ReeJl99Sf) . It is therefore important to carefully 
study these physical processes and determine to what extent 
each of them influences the characteristics of the host galax- 
ies. Nearby radio galaxies provide excellent opportunities to 
do this in detail. 

An important feature that is observed in a signifi- 
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cant fraction of active galaxies is the presence of out- 
flows of warm and hot gas from the central region (e.g. 
iHeckman et al.lll98ll : fVeilleux et all .2002.: iKrissl .20041 and 
references therein). In a number of nearby powerful ra- 
dio galaxies asymmetric emission-line profiles provide evi- 



10 3 km s" 1 (e.g. 


Tadhunterlll99ll:IVillar-Martin et all 


19991: 


Tadhuntcr ct al. 


2001 : iHolt. Tadhunter fc Morgantil 


120031: 


Tavlor, Tadhunter & Robinson 2003). Recently, a number 



of nearby radio galaxies have been found to contain such 
fast outflows not only in the ionized gas, but also in 
neutral hydrogen (H I) gas. Amo ng these radio galax- 
ies are 4C 1 2 50 (PKS 1345+12: IMorganti et al.1 l2004al) 
and 3C 305 iMorganti et alJ .20051) . An overview of the 
HI outflows found so far in radio gal axies is given by 
IMorganti. Oosterloo fc Tadhunterl J2004lh The exact driv- 
ing mechanisms for these fast outflows of neutral and 
ionized gas are not always clear, but they could con- 
sist of interactions between the ISM and the propagat- 
ing radio plasma (like the case of Seyfert g alaxy IC 5063; 
lOosterloo et all2000l IMorganti et al.l.2q04bl). AGN induced 
wind s (e.g. | Kr oll^!fc^ BegeliTuinr^86l Ir3alsara fc Krolikl 
Il993 : iDopita et alJ I2002T) or starbu rst-related phenomena 
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tral gas involved in the fast outflows puts serious constraints 
on the, most likely very energetic, driving mechanisms. 

An excellent case to study the physical mechanisms in 
detail is the nearby powerful radio galaxy 3C 293, in which 
we recently detecte d an outflow of neutra l hydrogen gas from 
the central region jMorganti et alj|2003l . hereafter Paper 1). 
Velocities of the outflowing HI gas reach ~ 1000 km s -1 , 
making it one of the most extreme cases of outflowing HI 
presently known. Unfortunately, due to limited spatial reso- 
lution of the observations we were not able to determine the 
exact location of this HI outflow. In this paper we analyse 
the characteristics of the ionized gas in 3C 293 and relate 
its kinematics with that of the outflowing neutral gas. We 
aim to gain a better understanding of the physical process 
that is responsible for the outflow and the effect it has on 
the ISM of the host galaxy. 

3C 293 is a powerful (Pi. 4 ghz ~ 2 x 10 25 W Hz" 1 ) 1 
edge- brightened Fanaroff-Riley type- II l|Fanaroff fc Rilevl 
1974) radio source with two ext ended radio lobes (e.g. 
iBridle. Fomalont fc Cornwelll Il98ll see also Figure . It 
also has a Steep Spectrum Core that consists of a 
true flat-spectrum radio core and a two-sided radio jet 
structure which is distorted and misaligned with re- 
spect to the outer radio lobes (e.g. lAkuior et alJ Il996t 
iBeswick. Pedlar fc Holowavl2002l see also Figure [TJ. 3C 293 
is hosted by the early type galaxy UGC 8782, that 
has a complex disk-like morphology with a bridge-tail 
structure that extends beyond a possible companion 
~ 33 arcsec to the southwest Jvan Breueel et alT Il984 
iHeckmanetalj Il986t lEvans et al.lll999l) " The spectroscopic 
work of Tadhun ter et al.l 820051) reveals post-starburst young 
stellar populations with ages between 0.1 and 2.5 Gyr 
throughout the galaxy. 3C 293 has a modest far-IR lumi- 
nosity of Z/fl r ~ 2.3 x 10 10 Lq. The galaxy contain s ex- 
tensive filamentary dust lanes (e.g. iMartel et al]ll99H) and 
large amounts of molecular and neutral hydrogen gas in the 
central few kpc (iBaan fc Haschick1ll98lt lH aschick fc Baanl 
[ 19851: lEvans et alJll999t IBeswick. Pedlar fc Holowavl l2002t 
IBeswick et alJl20 04). Given its peculiar optical morphology 
and the presence of large amounts of cold gas and young 
stars, 3C 293 has often been suggested to have been involved 
in a gas-rich galaxy-galaxy interaction or merger event, in 
which large amounts of gas have been funnelled into the nu- 
clear region and are now fuelling the AGN. The presence of 
a dense ISM, a young stellar population and a distorted ra- 
dio morphology make 3C 293 an excellent case to study the 
interplay between ISM, nuclear starburst and AGN activity. 



2 OBSERVATIONS 

Optical long-slit spectra were taken at the William Herschel 
Telescope (WHT 2 ) on May 12th 2001 (p. a. 60° and 135°) 
and in service mode on June 18th 2003 (p. a. 93°) using 
the ISIS long-slit spectrograph with the 6100A dichroic, the 

1 Ho = 71 km s -1 Mpc -1 used throughout this paper. At the 
redshift of 3C293, z = 0.045, this puts the galaxy at a distance 
of 190 Mpc and 1 arcsec = 0.92 kpc. 

2 The WHT is operated on the island of La Palma by the Isaac 
Newton Group in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 




Figure 1. Position of the slits along 60° (major axis host galaxy), 
135° (outer radio axis) and 93° (inner radio axis). The large- 
scale image is an overlay of WSRT radio contours (constructed 
from our 10 MHz data) superposed onto an optical DSS image 
of 3C 293. Contour levels in the WSRT image are 4.5, 6.5, 9.0, 
13, 18, 25, 35, 50, 70, 95, 130, 170 mjy beam" 1 . The zoom-in 
on the central region is an overlay of a MERLIN radio contin- 
uum image onto a HST F702W b roadband image (overlay from 
IBeswick. Pedlar fc Holowavl Eool see Section 13.11 for more de- 
tails). Contour levels in the MERLIN image are 9, 36, 144, 630 
mjy beam -1 . The regions C and E0 are the region of the core 
and the region of bright optical and radio emission 1 arcsec east 
of the core. These regions are described in detail in Section 13. II 



Table 1. Observational Parameters 







slit 


integration 






date 


p. a. 


width 


time (s) 


average 


seeing 




(°) 


(") 


blue red 


airmass 


(") 


5/12/2001 


60 


1.3 


2400 2400 


1.085 


~ 0.8 


5/12/2001 


135 


1.3 


1200 1200 


1.037 


~ 0.8 


6/18/2003 


93 


1.03 


3600 2400 


1.02 


~ 1.0 



R300B and R316R gratings at the blue and red arm and the 
GG495 blocking filter on the red arm to cut out second order 
blue light. This resulted in a wavelength coverage from about 
3500 to 8000 A. The slit was aligned by eye and centred on 
the brightest part of the galaxy. The different slit positions 
are indicated in Figure Q Table Q gives a summary of the 
observational parameters. 

We used the Image Reduction and Analysis Facility 
(IRAF) to reduce the data. A description of the data re- 
duction of the May 1 2th 2001 data is already given by 
iTadhunter et alJ J2005t> . For the June 18th 2003 data (p. a. 
93°) the standard calibration (bias subtraction, flatfielding 
and wavelength calibration) was done. We used dome and 
arc exposures taken at approximately the same position and 
time as the source spectra. We took out a significant tilt 
of the slit in the spatial direction by applying a 1st order 
correction using two stars that were in the slit. With the 
tilt removed, the spectra are aligned within one pixel. Af- 
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ter the background subtraction the frames were combined 
and cosmic rays removed. Two standard stars (BD+26 2606 
and BD+33 2642) were used for the flux calibration. The 
wavelength calibration was checked for both the 2003 and 
2001 data using night skylines in spectra that were reduced 
without background subtraction. We corrected for any sys- 
tematic off-set of the lines in the different spectra (likely 
due to flexure in the spectrograph) during the analysis of 
the spectra. The resulting accuracy of the wavelength cali- 
bration is within 0.1 A for all the red spectra and 0.4 A for 
all the blue spectra, with no obvious systematic errors. The 
A-resolution of the spectra is ~ 4A. 

For the analysis of the spectra, we used the Starlink 
package FIGARO to bin consecutive rows of pixels into a 
one-dimensional spectrum. In this way, we produced a series 
of spectra across the slit. These spectra were further anal- 
ysed with the Starlink package DIPSO by fitting Gaussian 
profiles to the spectral lines. Velocities used in this paper 
are heliocentric velocities. 

In order to study neutral hydrogen gas in emission 
around 3C 293 we obtained Westerbork Synthesis Radio 
Telescope (WSRT 1 ) observations on May 7th 2001, using the 
10 MHz band. The goal of this study was to look for pos- 
sible HI structures that could be reminiscent of a merger 
event in 3C 293. In the line data we detect three nearby 
companions of 3C 293 in HI em ission, which we describe 
in detail inlEmo nts et all <l2004h . However, as explained in 
lEmonts et al.l i2^0^ ~the bandwidth of these observations 
turned out to be not ideal for looking for H I emission di- 
rectly associated with the host galaxy of 3C 293. In the 
present paper we use the WSRT data to make a new contin- 
uum image of 3C 293. The channels containing any deep H I 
absorption were excluded from the continuum fitting. We 
used the MIRIAD software to make the continuum image 
shown in Figure Q which we use as a reference to show the 
slit alignments. The detection level of the radio continuum 
is limited by a dynamic range ~ 500, while the peak inten- 
sity is 3.6 Jy beam" 1 . The beam-size is 25.3 x 11.9 arcsec at 
p.a. 1.1°. The total power is Pi. 4 gh z w 2 X 10 25 W Hz" 1 . 



3 RESULTS 

Extended continuum and emission-lines are detected in all 
our optical spectra of 3C 293. Figure|5|shows a typical spec- 
trum of 3C 293 in the regions E0 and C (these regions are 
indicated in Figure Q and will be described in detail below). 
In this Section various features of the emission-line gas in 
3C 293 along the different slits are presented. 

For the analysis of the kinematics of the gas we will 
mainly use the [S II] and [O II] lines. These are among the 
stronger emission lines detected and they do not suffer from 
absorption features in the underlying continuum as can be 
the case for the stronger Ha line. 



The WSRT is operated by the Netherlands Foundation for Re- 
search in Astronomy with support from the Netherlands Founda- 
tion for Scientific Research. 
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Figure 2. Spectrum of the regions E0 and C. The spectrum is a 
composite spectrum of both the red and the blue spectrum along 
p.a. 93° . The flux is summed over an aperture of 0.8 arcsec. See 
Section l.3.3l for more details on these spectra. 



3.1 Gas kinematics along the inner radio axis 

In this section we concentrate on the analysis of the kine- 
matics of the gas along the inner radio axis, i.e. p.a. 93°. 
Along this position angle the ionized gas clearly shows the 
most complex kinematics. This is illustrated by the 2D spec- 
trum of the Ha+[N II]AA6548,6583 region shown in Figure 
|3] In this figure, the locations of the various features de- 
tected in the optical spectrum are compared with features 
seen in the radio and the HST image. Figure |H] (bottom) 
shows a MERLIN radio contour map overlaid on an opti- 
cal HST image of the inner 5 kp c of 3C 293 (figure from 
iBeswick. Pedlar fc Holowavl 120021 also shown in Figure 0. 
The accuracy of the overlay is within 0.3 arcsec. The HST 
WFPC2 image was taken with a broadband F702W filter 
and includes continuum in the wavelength range just below 
the strong Ha+[NII] lines. In Figure (middle) an inte- 
grated brightness profile from the HST image is taken across 
the region that is covered by the 93° slit, as indicated by the 
dotted region in the HST image. Figure [5] (top) shows that 
the different features in the HST image can be identified in 
the red part of our 2D-spectrum along p.a. 93°. Several re- 
gions of interest are indicate d. Region C, when ob served in 
the radio at VLBI resolution iBeswick et aljEo04h . includes 
the actual core as well as a radio knot about 0.4 arcsec west 
of the nucleus. Region E0 is the region of peak intensity in 
the optical, about 1 arcsec east of the nucleus. This is also 
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Figure 3. Bottom: MERLIN radio contour map overlaid on an 
optical HST image of the inner f ew kpc of 3C 293 (figure from 
iBeswick. Pedlar fc Holowavll20o3l Contour levels are 9, 36, 144, 
630 mjy beam . Middle: Integrated profile from the HST im- 
age across the region that is covered by the 93° slit (indicated 
by the dotted region in the bottom plot). Top: Spectrum along 
p.a. 93° (Ha+[NII]). Contour levels: 0.75, 1.67, 3.4, 6.1 X 10~ 17 
erg cm - 2 s — 1 A — 1 . The features in the spectrum can be traced in 
the HST image. 



the region in which the radio continuum peaks. Regions El 
and Wl are regions east and west of the nucleus just at the 
edge of where fainter radio continuum is detected. Further 
out to the east and west (in regions where so far no radio 
continuum has been detected) are the regions E2 and W2. 
In Figure 2] we show the 2D-spectrum of the [S II] dou- 



blet. Also shown in Figure 2] are fits to the [S II] line at 
various places along the slit. As mention above, the [S II] 
line is used for the kinematic analysis in order to avoid pos- 
sible effects of absorption features that may affect the Ha 
line. Moreover, the region of the nucleus is more clearly seen 
in the red part of the spectrum than in the blue part (due 
to the higher sensitivity of the ISIS spectrograph in the red 
as well as complex extinction ac ross the spectral ra nge due 
to prominent central dust-lanes: iMartel et al.lll999T) . which 
favours the use of [S II] above [O II] for the detailed kinematic 
analysis presented in this Section. 

In Figure |1] we can immediately identify a narrow com- 
ponent in most of the fits. To fit the narrow component of 
the [S II] doublet we use two Gaussians with equal width and 
with wavelength separation constrained, following the phys- 
ic al parameter s (corr ected for the redshift of 3C 293) given 
bv lOsterbrockl il989l) . At most places along the slit we need 
a second, broader component to get a good fit (again this is 
a doublet component with equal width and wavelength sepa- 
ration constrained). Figure [S] shows various results from the 
line-fitting of both components of the [S II] line in the spec- 
trum along 93° . Plotted against the distance of the features 
from the nucleus are the velocities of both the narrow and 
broad component relative to the systemic velocity of 3C 293, 
the velocity shift between the narrow and broad component, 
and the width of the broad component along the slit. Also 
shown is the line flux of both the narrow and broad com- 
ponents. Table |21 summarizes the kinematic information of 
both components in the different regions. 

From Figures [I] and |H| and Table [5] we can distinguish 
several features with distinct gas kinematics. In Figurc|5Ja), 
the narrow components (connected with the dashed line) 
trace what appears to be a regular rotating feature that we 
identify with the l arge s cale gas disk already observed by 
Ivan Breueel et all il984T) . We will discuss more about this 
disk in Section 13.21 In the following we briefly discuss the 
characteristics of the ionized gas in each region: 

• Region C: Region C is the region that includes the 
nucleus. In region C we need to include a broad component 
to get a good fit to the emission lines (Figure This broad 
component, however, is strong compared to the broad com- 
ponent at other places along the slit and it peaks at the 
same velocity as the narrow component. It therefore might 
reflect a more turbulent state of the ionized gas close to the 
nucleus of the galaxy or the presence of a circum-nuclear 
disk. Similar phenomena have been suggested to to explain 
a broadening of t he deep H I absorption profile in the nuclear 
region of 3C 293 jBeswick. Pedlar fc Holowavll2002l) . 

At the location of the nucleus, we very accurately deter- 
mine the systemic velocity of 3C 293 to be v sya = 13, 450T35 
km s" 1 or z = 0.04486 ± 0.00012, which is the redshift of 
the narrow component in region C. The error has been de- 
termined from the uncertainty in fitting the Gaussian com- 
ponents, from the inaccuracy of the wavelength calibration 
and from the uncertainty in the exact position of the nucleus 
in the optical spectrum. 

• Region EO: Region E0 is the region about 1 kpc east 
of the nucleus, where both the optical as well as the radio 
continuum peak in intensity. In region E0 we detect, in addi- 
tion to a narrow component, a broad component of the [S II] 
line that is highly blueshifted with respect to the velocity 
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Figure 4. Spectrum of the [S II] AA6716,6731 line (middle; contour levels 0.2, 0.48, 0.72, 1.1, 1.6, 2.5, 3.6, 4.6 X 10~ 17 erg cm^s -1 A -1 ). 
Also shown are the fits to the [S II] line at various places along the slit. At most places along the slit we need two components to get a 
good fit. The apertures used to make the ID-spectra that we present here are written in the bottom-right corner of the plots (note that 
the plot of region W2 is taken with an aperture of 0.8 arcsec, but a similar spectrum is seen throughout the region W2 as it is indicated 
in Figure 15*1 , 



Table 2. Kinematic information of the various narrow and broad components observed in the different regions. Given is the FWHM 
(corrected for instrumental broadening) of the fit to both components as well as the velocity shift (Av) between the centre of the broad 
component and the centre of the narrow component in each region (negative shift means that the broad component is blueshifted). Also 
given is the [S II] A6716+A6731 flux of the various components (given as the calculated mean flux over an aperture of 2 arcsec). Due to 
the faintness of the narrow component across region W2 we could not obtain a single good value for the width of this narrow component, 
therefore we left it out of this Table. JVote that the broad component represents different gas features in the different regions (see text). 

Region FWHM„ arr ™ FWHM 6roati Av Flux narrotu Flux 6rotld 

(km s" 1 ) (km s -1 ) (km s" 1 ) (xlO -16 erg s _1 cm -2 ) (xlO -16 erg s _1 cm -2 ) 

E2 - 664 ± 27 -202 ± 22 - 1.00 ± 0.04 

El 248 ± 15 687 ± 80 -114 ± 48 0.90 ± 0.10 0.61 ± 0.12 

E0 341 ± 10 729 ± 127 -542 ± 56 4.02 ± 0.20 1.35 ± 0.26 

C 411 ± 20 1032± 57 26 ± 31 3.12 ± 0.39 5.73 ± 0.64 

Wl 112 ± 55 398 ± 68 115 ± 75 0.19 ± 0.12 0.30 ± 0.12 

W2 - 264 ± 40 152 ± 33 - 0.41 ± 0.09 
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Figure 5. Line-fitting results of the [S II] line in the spectrum along 93°. In the various plots the small squares (connected with a dotted 
line) represent the narrow component of the gas. The circles represent the broad component. Vertical error bars represent uncertainties 
in the fits and in the A-calibration; horizontal error bars represent the width of the used apertures. The horizontal axis is centred on 
the nucleus. Plotted against the spatial off-set from the nucleus are: (a) the velocity of both the narrow and the broad component w.r.t. 
the systemic velocity (v syB = 13, 450 ± 35 km s — x ; see Section l3.ll : (b) the velocity shift between the broad and narrow component (the 
dotted line, velocity shift = 0, corresponds to the velocity of the narrow components); (c) the width of the broad component; (d) the 
[S II] A6716+A6731 line flux (mean value over an aperture of 0.2 arcsec) of both the narrow and broad component. 



of the narrow component and the systemic velocity of the 
system. A similar component is seen in the [O II] lines. The 
FWHM (full-width at half the maximum intensity) of this 
component is 729 ± 127 km s _1 (corrected for instrumental 
broadening) and its centre is blueshifted by 542 ±56 km s _1 
w.r.t. the velocity of the narrow component that traces the 
regular rotating gas disk (see Figure^. The width of the 
broad component and the amplitude of the velocity shift 
(w.r.t the narrow component) are so large that they cannot 
be due to gravitational motions of the gas. Thus, we inter- 
pret this broad, blueshifted component as a fast outflow of 
ionized gas. The exact extent of this outflow - both along 
the radio jet as well as perpendicular to it - will be discussed 
in Section f4.ll 

• Region El: Region El is the region at the edge of the 
weak radio emission about 4 kpc east of the nucleus. In re- 
gion El the [SII] emission line shows a broad component 
with FWHM = 708 ± 80 km s" 1 . In this case, the centre of 



this broad component is blueshifted by only 114 ±48 km s - 
w.r.t. the narrow component. Thus, also this component 
may represent a gas outflow, although not as pronounced 
as in region E0. 

• Region E2: A very remarkable feature can be seen in 
region E2 at ~ 7 arcsec east of the nucleus. In Figure |!J] this 
feature clearly has a "comma" -like shape and the ionized 
gas disk seems distorted. This is strengthened by the fact 
that in region E2 we cannot identify the narrow component 
in the [S II] line. Instead, we fit a single broad component in 
region E2. The FWHM of this feature is 664 ± 27 km s" 1 . 
When extrapolating the velocity that the disk would have 
in region E2 (as traced by the narrow component in the 
regions adjacent to E2; see Figure |HjA), the centre of this 
broad component is blueshifted by 202 ± 22 km s _1 , again 
indicating the presence of gas outflow also in this region. In 
region E2 there is no detection of radio continuum down to 
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Figure 6. Position-velocity diagram of the [Oil] line along the 
major axis of the galaxy (p. a. 60°). The squares trace the narrow 
component; the circles the broad component. Point (0,0) repre- 
sents the nucleus. The off-set on the horizontal axis is the off-set 
(in arcsec) from the nucleus (error ± 0.6 arcsec). On the vertical 
axis the velocity is plotted w.r.t. the systemic velocity that we de- 
rive at the nucleus (v S j /s =13,450 ± 35 km s ; see Section 13.11 . 
The vertical error bars represent uncertainties in the Gaussian 
fits and inaccuracies in the wavelength calibration. 



a limit of 3 mjy beam 1 in radio continuum maps presented 
by e.g. lBridle. Fomalont fc CornwelJ l|l98l|) . 

• Regions Wl & W2: Regions Wl and W2 are the 
regions west of the nucleus. In regions Wl and W2 the 
emission-line profile, although very faint, is clearly asym- 
metric. By including an additional, broader component we 
get an excellent fit to the emission-line profile in this part 
of the galaxy. The narrow component on the western side 
nicely traces the velocity curve of the gas disk, while the 
broader component is redshifted with respect to this narrow 
component. However, the FWHM of the fit to this broader 
component in regions Wl and W2 is only 398 ± 68 (Wl) 
and 264 ± 40 (W2) km s _1 . Furthermore, the shift between 
the narrow and broad component is only 115 ±75 (Wl) and 
152 ± 33 (W2) km s . This redshifted component is not as 
pronounced as the blueshifted component seen on the east- 
ern side of the galaxy. It might, therefore, indicate that only 
a mild outflow is present on the western side of the galaxy 
(with gas moving away from us). Another possibility is that 
this component might be caused by line-of-sight effects when 
we look through the almost edge-on disk of the galaxy. The 
asymmetric emission-line profile would then represent the 
narrow component gas at various locations in the disk, in- 
tercepted by our line-of-sight. With the data available to us 
we cannot disentangle these two possibilities. 

The main result from this Erst part of the kinematic 
study of the ionized gas is that there are a number of loca- 
tions along the position angle of the inner radio jet where gas 
outflows are observed. In addition to this, the most promi- 
nent outflow is observed at the location of the stronger radio 
continuum (E0). This will be further discussed in Section^ 



3.2 Extended gas disk 

The spectrum taken along p. a. 60° follows the major axis of 
the galaxy. Therefore we also analysed the kinematics of the 
emission-line gas in this direction. Figure|5]show a velocity- 
curve of the ionized gas, constructed by fitting Gaussian 
profiles to the [O II] line, using apertures of 0.6 arcsec at 
various places along the 60° slit. Being one of the stronger 
lines, [O II] can be traced out to a larger distance from the 
centre than other lines. From this diagram it is clear that the 
[O II] emission-line gas fo llows the kinematics of a regular ro- 
tating disk (also seen bv lvan Breueel et al.l [l984 % ) . The plot 
also clearly shows how the blueshifted component observed 
in region E0 (the only one included in this slit position) is 
well separated from the gravitational motion of the regu- 
lar rotating disk gas. The total diameter of the gas disk is 
30 kpc, which is larg e r tha n the 19 kpc disk already found 
bv lvan Breueel et alJ dl984l) . A small difference in slit align- 
ment might exp lain the slight deviation o f our velocity curve 
from the one of Ivan Breueel et alJ lll984l) . 

We centre the velocity curve on the systemic velocity 
that we derived in Section 13. II («sys = 13, 450 ± 35 km s _1 ). 
Our value for v sys is in agreement with v sya = 13, 478 
km s -1 , fou n d wi th HI absorption at the core by 
iBeswick et alJ i2004l) . The redshift velocity of the disk gas 
in region E0 is roughly 13,500 km s _1 , consistent with 
the velocity of the deep, narrow H I absorption that is 
located against the enhance d radio emission in reg i on E0 
iBeswick. Pedlar fc Holowavl 12001 IBeswick et alJ 12004) 
and with CO measurements T llEvan^^uTTr999l) . From 
our rotation curve along p. a. 60° we derive a velocity 
gradient of ~ 57 km s _1 arcsec -1 in the inner region of 
the galaxy. This i s cons istent with the gradient seen by 
Ivan Breueel et ail {l984) and with the gradient of ~ 46 
km s -1 arcsec -1 derived f r om H I absorption studies by 
IBeswick. Pedlar fc Holo'rovl <2002l) . 



3.3 Ionization and density of the gas 

Table gives the values of the fluxes of the narrow compo- 
nent of various emission lines measured in the regions E0 and 
C in the spectrum along p. a. 93° (see Figure The fluxes 
are the total emission-line fluxes, estimated from Gaussian 
fits, across an aperture of 0.8 arcsec (corresponding to a 
spatial binning of 4 consecutive pixels in the 2D-spectra) 
centred on regions E0 and C. 

The line-ratios in Table 3 are not corrected for extinc- 
tion. The ratio Ha/H/3 = 9.5 in region E0 is about a fac- 
tor 3 too high when comparin g it with the hyd rogen re- 
combination analysis done by lOsterbrockl J1989T) . There- 
fore it is clear that H/3 is strongly affected by redden- 
ing (caused by variou s dust lanes that cr oss the inner re- 
gion of 3C 293, e.g. iMartel et aT] Il999l) and/or by un- 
derlying H/3 absorp tion associated with t he young stel- 
lar population that iTadhunter et all {2005) found in this 
galaxy. Furthermore, the Ha emission line could be af- 
fected by underlying atmospheric absorption. All this makes 
the analysis of the line ratios more complicated. Never- 
theless, when comparing with diagnostic diagrams (e.g. 
[ Baldwin. Phillips fc Terlevichl Il98ll : IVeilleux fc OsterbrocM 
Il987t iDopita fc Sutherland! Il995f) . we can use some of the 
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Table 3. Fluxes of the narrow component emission-line gas - rel- 
ative to H/3 - across an 0.8 arcsec aperture in regions E0 and C. 
The H/3 flux is given in erg s -1 cm -2 . The errors are also scaled 
to the value of the H/3 flux and reflect both the uncertainty of the 
Gaussian fits as well as the inaccuracy in determining the contin- 
uum level. Also included are the fluxes of the broad component 
of the [O II] and [S II] gas (scaled to the value of the narrow com- 
ponent) in region E0. 





Region EO 


Region C 


Element 


Flux 


Flux 


(narrow comp.) 


relative to H/3 


relative to H/3 


H/3 flux 


2.32±0.24xl(r 16 


1.15±0.10xl0- 16 


[OII]A3727 


4.2 ±0.2 


4.9 ± 0.1 


[OIII]A4363 


<0.3 




[Offl]A5007 


1.2 ± 0.2 


2.6 ± 0.1 


[OIJA6300 


1.2 ±0.4 


5.2 ± 0.1 


H Q 


9.5 ±0.4 


15.7 ± 0.7 


[N II]A6583 


7.7 ±2.9 


20.4 ± 3.8 


[SII]A6716±A6731 


7.4 ±0.3 


21.9 ± 0.3 


[O II] broad/ narrow 


0.32 ±0.05 




[S II] broad /narrow 


0.30 ±0.05 





lines, which lie close together so that the line ratios are not 
affected by extinction too much (e.g. [Nll]/Ha, [SII]/Hq 
and [Ol]/Ha ratios, together with [O III] /H/3), as good di- 
agnostics for the classification. Using these diagnostics we 
can conclude that region E0 in 3C 293 has a LINER-type 
Low Ionization Nuclear Emission-line Region) spectrum 
Heckman"lll98Cf) . Thus, the narrow component gas in re- 
gion E0 is likely either shock ionized or photo-ionized by 
the AGN at low ionization parameter. The line ratios are 
inconsistent with photo-ionization by O and B stars (HII 
region) or by a strong power-law AGN component at high 
ionization parameter. 

The [O II] and [S II] fluxes of the broad component gas 
associated with the outflow in region E0 are also listed in 
Table [3] We were not able to make an accurate fit of the 
broad component in region E0 to any other line than the 
strong [O II] and [S II] lines. Other lines also show evidence 
for a broad component, but due to the underlying contin- 
uum and the limited signal-to-noise (in combination with the 
complexity of some of the emission lines, like e.g. [Nll]±Ha) 
we were not able to get a reliable fit to these broad compo- 
nents. 

The emission-line gas in region C (fitted with a single 
Gaussian component) appears to have ionization character- 
istics similar to the narrow line gas in region E0, although 
the measured Ha/H/3 = 15.7 indicates that region C is likely 
even more affected by dust-extinction than region E0. In re- 
gion C the gas is also likely ionized either by shocks or by 
the AGN at low ionization parameter. 

Electron densities o f the gas can be derived from the 
[SII]A6716/A6731 ratio llOsterbrocklll989l) . Figure □ shows 
[S II]A6716/A6731 for both the narrow and the broad com- 
ponent in the various regions along p. a. 93° (the fit is too 
uncertain across region W2 and therefore we leave it out 
of the plot). The broad component represents different fea- 
tures in the various regions (Section 13. II . which might be 
physically unrelated. Nevertheless, given the relatively large 
uncertainties, we do not see any difference in the [S II] ratio 
in the various regions. The weighted mean [S II]A6716/A6731 
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Figure 7. [S II]A6716/A6731 of both the narrow emission-line 
gas component (closed squares) as well as the broad emission- 
line gas component (open circles) in the various regions along 
p. a. 93°. Shown are the dotted lines that correspond to densities 
of 10 and 300 cm" 3 (assuming T = 10,000 K, see lOsterbrockl 
Il98fll) as well as the maximum and minimum value for the [S II] 
doublet line ratio that is physically possible (the upper dashed 
line corresponds to n c < 10 1 cm -3 , the lower dashed line to 
n c > 10 5 cm -3 . The error bars reflect the uncertainty in fitting 
the Gaussian components. 

ratio in Figure[7|is 1.23T0.08 for the narrow component and 
1.19 ± 0.07 for the broad component (excluding the broad 
component in region C, that does not represent an out- 
flow; see Section ^. H . This corresponds to densities of about 
2 ± 1 x 10 2 cm -3 for the narrow component and 3 ± 1 x 10 2 
cm -3 for the broad component gas, assuming T = 10, 000 K 
(since [O III]A4363 is not detected we are not able to derive 
the electron temperature of the emission- line gas). Thus, we 
do not see any significant difference between the density de- 
rived for the narrow and the broad component gas. 



4 FAST GAS OUTFLOW IN REGION EO 

The fastest outflow of gas is seen in region E0. In this region 
the emission lines are particularly strong and they allow 
some further analysis of the characteristics of this outflow. 
In particular we investigate in this Section the spatial extent 
of the region showing the fast outflow and we estimate the 
mass of the outflowing gas. This information will be used 
later on for building a possible scenario about the origin of 
this outflow. 

4.1 Spatial extent of the gas outflow in region EO 

In Figure |8] we analyse the extent of the region in which we 
detect the fast outflow. We do this by tracing the broad, 
blueshifted component of the [Oil] emission line along the 
different position angles (since the [O II] line is stronger 
than the [S II] line it provides a better indication for the 
presence of the broad component). The top plots show in 
contours the HST F702W image of the central few kpc of 
3C 293. The radio continuum of the MERLIN map (from 
iBeswick. Pedlar fc Holowavl [2002) is plotted in grey-scale 
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Figure 8. In the top plots contours of the HST F702W image of the inner region of 3C 293 are overplotted onto a grey-scale plot of 
the MERLIN radio continuum map. The dotted line presents the rotation axis of 3C 293 through the nucleus. The boxes indicate the 
region in which the broad, blueshifted component of the [O II] emission line is detected (see text for the accuracy of the position of the 
boxes). In the plots below that, various observed parameters are plotted. The first plot shows the width of the broad (circles) and the 
narrow (squares) component (corrected for instrumental broadening). The plots below that show the velocity shift between the broad 
and narrow component. The bottom plots show the line flux of both the broad (circles) and narrow (squares) component. The off-set on 
the horizontal axis is w.r.t. the nucleus, which we can determine with an uncertainty of 0.6 arcsec for the [Oil] line. The horizontal bars 
in the bottom-right corner of the plots indicate the apertures over which we summed rows of pixels to extract the ID-spectra (note that 
for the line flux along p.a. 60° and 135° we used a 0.4" aperture and a 1.3" slit width, while for the p.a. 93° the aperture was only 0.2" 
and slit width 1.03"). 



© 2004 RAS, MNRAS 000,0-?? 



10 B.H.C.Emonts et al. 



and clearly shows the knots of enhanced radio emission in 
the regions C and EO. The length of the boxes indicates 
the region in which the broad, blueshifted component is de- 
tected; the width of the boxes corresponds to the width of 
the slit. Note that the slit was aligned by eye to cover the 
brightest region and Figure 03 merely shows the most likely 
positioning of the slit. The uncertainty in the position of the 
nucleus in the spectra around the [Oil] line (the "0" -point 
for the spatial off-set) is 0.6 arcsec. Also shown in Figure 
[H] are the width of the broad and narrow component, the 
velocity shift between the two components and the line flux 
of both components across the region in which the broad 
component is detected. 

From Figure|H]and the seeing listed in TableQit is clear 
that the region in which the broad, blueshifted component 
is observed, i.e. the fast-outflow region, is spatially resolved. 
Not only is it resolved along the direction of the inner radio 
axis (where, moreover, we know already other regions of gas 
outflow do exist; see Section f3.H . but also perpendicular to 
it. The gas outflow is observed in a region of 2 x 1.5 kpc 2 . 



4.2 Mass of the gas outflow in region EO 

The mass of the gas in a line emitting region is related to 
the H/3 luminosity in that region by: 



M eas = m D 



L(R/3) 



(1) 



(see IOsterbrocklll989D . where n c is the electron density 
(cm -3 ), m p the mass of a proton (kg), L(H/3) the H/3 lu- 
minosity (erg s _1 ), the effective recombination coeffi- 
cient for H/3 (cm 3 s _1 ) and hvnp the energy of a H/3 photon 
(erg). Assuming that the fast outflow is spread over an area 
of 2 x 1.5 kpc 2 (Section 14.111 and accordingly over 1.5 kpc 
in our line-of-sight direction, we first estimate the luminos- 
ity of the narrow component gas of the disk in this outflow 
region to be -L(H/3) narl . ow = 1.3 x 10 40 erg s _1 . Here we as- 
sume that the extinction corrected flux (about 3x the flux 
given m Table 03 which was measured with a slit aperture 
of 0.8 x 1.03 arcsec 2 ; see Section f3. 311 . uniformly covers the 
region of outflow. This is a reasonable assumption, as can 
be seen from Figure |S| (bottom). To estimate the mass of 
the outflowing gas we need L(H/3) of the broad component 
gas. Unfortunately, the quality of the data did not allow us 
to unambiguously fit a broad component to the H/3 line (see 
Section ^. 31 . Nevertheless, when assuming that the line ratio 
HPbroad/Hflnarrow ~ 0.3 (similar to what we measured for 
[Oil] and [S II] ; Table 0U, we derive for the H/3 luminosity 
of the broad component L(H/3)broad ~ 4.3 x 10 39 erg s _1 . 
Using the density derived in Section 13.31 we estimate that 
the total mass of the outflowing ionized gas in region EO is 
M » 1 x 10 5 Mq . The total ionized gas mass for the narrow 
component disk gas in the outflow region is M w 6 x 10 5 
M e 

With the oversimplified assumption that the gas is 
locked up in clouds with a derived density of 3 x 10 2 cm -3 
one can make an estimate for the filling factor (Ff) of these 
clouds over the region where we detect the fast outflow: 



where M is the total mass of the outflowing gas, n the den- 
sity of the gas, m p the mass of a proton and V the vol- 
ume of the region where the fast outflow is detected. Us- 
ing the derived properties of the gas involved in the fast 
outflow, Equation [5] gives a filling factor of the outflow- 
ing gas of Ff ~ 3 x 10 -6 . Of course this is a rough es- 
timate, since the distribution of the gas is unknown, but 
it is consiste nt with values derived in other powerful ra- 
dio galaxies Ivan Breueel et alJll985l Il986t Uleckman et alJ 
119891: iKoekemoer fc Bicknefll ll99Sft . It shows that we are 
dealing with a clumpy medium in the central region of 
3C 293. 



5 COMPARISON WITH THE H I OUTFLOW 

In Paper 1 we presented evidence for a fast outflow of neutral 
hydrogen gas in 3C 293. Figure [5] (from Paper 1) shows the 
broad but shallow HI absorption, mainly blueshifted w.r.t. 
the systemic velocity, that represents this fast outflow of 
HI. This absorption is detected against the central radio 
continuum using the 20 MHz band at the WSRT, but due 
to the low spatial resolution of the WSRT we are not able 
to locate the region in the inner radio structure from which 
the HI outflow originates. 

In lEmonts et al] (120041) we made a Gaussian fit to this 
broad HI absorption and derived a FWHM of 852 ± 41 
km s _1 . The centre of this Gaussian fit is 345 ± 24 km s _1 
blueshifted from the deepest, narrow HI absorption feature, 
which traces quiescent gas like l y located in an extended H I 
disk llHaschick fc BaaiJ Il985| : iBeswick. Pedlar fc Holowavl 
2002t lBeswick et al B200 4) . The characteristics of the broad, 
blueshifted H I component look similar to those of the broad 
component of the ionized gas observed in region E0. A com- 
parison between the two is shown in Figure 1101 (from Pa- 
per 1). Region E0 therefore appears to be the most likely 
location for the fast HI outflow. The fact that the radio 
continuum is strongest in region E0 further strengthens this 
conclusion. The observed outflow of HI cannot be associ- 
ated with the outflow of ionized gas that we see in region 
E2, since no strong radio continuum is detected in region 
E2. Moreover, at the resolution of the WSRT observations, 
the outflow would appear not coincident with the nucleus 
but displaced with respect to it. 

The similarity between the characteristics of the HI out- 
flow and the ionized gas outflow in region E0 suggests that 
they likely originate from the same driving mechanism. 

If the broad component of the HI absorption is indeed 
located against the region E0, we can derive a new value for 
the H I column density. The peak intensity of the bright ra- 
dio emission in region E0 in the MERLIN map of FigurcQis 
1.29 Jy. Therefore, the optical depth that the broad HI ab- 
sorption feature has if it is located in region E0 is r = 0.38%. 
This corresponds to column densities of the outflowing HI 
gas of ~ 6 x 10 20 cm -2 assuming T sp i n = 100 K. This value 
for the column density likely represents a lower limit, as in 
the extreme conditions under which these fast outflows oc- 
cur the T SD i n is likely (much) higher (e.g. iBahcall fc Ekersl 

36). 



Ff 



M 



nm p V 



(2) 



ll969tlMalonev. Hollenbach fc Tielenslll99i 

Assuming that the outflowing HI gas extends over the 
same region as the fast outflow of ionized gas, the H I column 
in front of the radio plasma will extend over ~ 0.75 kpc. 
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Figure 9. Right: Continuum image of 3C 293 made with our 10 MHz data. Left: Zoom-in on the HI absorption spectrum detected 
against the unresolved central radio continuum using the 20 Mhz band (Figure from Paper 1). Indicated is the full width at zero-intensity 
of the H I outflow. 




Figure 10. [Oil] emission line in region E0 and HI absorption 
profile (inverted) plotted in one image. 



The mean density of the neutral gas involved in the fast 
outflow will therefore be ~ 0.3 c m" 3 , which is comparable 
to the densities of ~ 0.1 cm" 3 that lHaschick fc Baanl ijlflSflT l 
found for HI clouds that are falling onto the AGN. However, 
given the very low filling factor that we derive for the ionized 
gas (Section I4.2|l . the neutral gas will most likely also be 
locked up in a clumpy medium with locally much higher 
densities. If the outflowing HI gas extends over the same 
region as the fast outflow of ionized gas in region E0 we 
derive a total HI mass of A/hi ~ 10 7 A/0 involved in the 
fast outflow. Although uncertainties are large (mainly due to 
the uncertainty in the exact location and extent of the HI 
outflow and due to the unknown spin temperature of the 
gas), the bulk of the outflowing gas nevertheless appears 
to be in a cold, neutral state and only a minor fraction is 
ionized. 



6 DISCUSSION 

The long slit spectra presented in this paper show complex 
kinematics of the ionized gas in the radio galaxy 3C 293, in 
particular along the position angle of the inner radio jets. 
At a number of locations we find outflows of ionized gas. 
The broader and most blueshifted gas-component, that we 
identify with a fast (~ 1000 km s" 1 ) outflow of ionized gas, 
is detected in the region E0, close to the brighter radio emis- 
sion found along the jet. This is intriguing and suggests that 
the presence of the radio jet may play a role in producing 
the outflow. The total mass of warm gas involved in the fast 
outflow is ~ 1 x 10 Mq. There are more regions in which 
complex kinematics of the ionized gas are seen, although less 
extreme than what we detect in region E0. 

In addition to this, one of the most interesting results 
obtained is the similarity in the kinematic properties be- 
tween the ionized gas outflow in the region E0 and the out- 
flow of neutral hydrogen presented in Paper 1. This suggest 
that the two outflows are coming from the same region (that 
would therefore be about 1 kpc from the nucleus) and that 
they are likely driven by the same mechanism. The total H I 
mass of 10 7 Mq involved in the outflow shows that most 
of the outflowing gas is in its neutral state and that only a 
small fraction is ionized. 

In the following, we will discuss some of the possible 
mechanisms that can cause such (off-nuclear) outflows, fo- 
cusing in particular on the case of the region E0, for which 
detailed information could be derived from the optical data. 
Although the scenario of a jet-ISM interaction appears to be 
the most likely mechanism for the fast outflow in region E0, 
other mechanisms are also possible. We also briefly discuss 
what could be the cause for the milder outflows observed in 
other regions across 3C 293. 
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6.1 Driving mechanism for the fast gas outflow in 
region EO 

Here we consider three possible mechanisms to explain the 
detected gas outflows, in particular the fast outflow detected 
in region EO. 

• Starburst driven wind: 

A possible driving mechanism for the outflow is ac- 
celeration of the gas by winds from young stars. By 
fitting the Spectral E nergy Distribution (SE P) at vari- 
ous places in 3C 293, iTadhunter et al] (|2005|) detected a 
young/intermediate age stellar population of 0.1 - 2.5 Gyr 
throughout 3C 293. Although starburst winds can gen- 
erally produce outflows up to 1000 km s _1 , can easily 
carry enough energy for an outflow such as observed in 
3C 293, and can predominantly shock-ionize the gas (e.g. 
iHeckman. Armus fc Milevlll990ft . it is not clear whether a 
starburst wind can survive over a period of 0.1 - 2.5 Gyr. 
Even if this is possible, a wind from this young/intermediate 
age stellar population would have to be "fossil" in order to 
still be visible at present. An outflow driven by such a fossil 
wind would no longer be localized to the region E0 where 
we see it, but would be many tens of kpc in extent. 

It should be noted, however, that bright complexes 
of UV emi ssion have been found throughout the disk 
of 3C 293 jAllen et al] I2002T) . Some of these complexes 
could contain young O and B stars, in which case winds 
from these young stars may contribute in some regions to 
the complex kinematics of the gas that we observe. The 
relatively poor accuracy of the astrometry between the 
UV image and the radio data does not allow us to reliably 
estimate how these complexes are located with respect to 
the radio features. It is worth noting, however, that one 
of the brighter UV complexes (one with an intriguing V- 
shaped morphology) could be located close to the region E0. 

• AG N radiation: 

iKrolik fc Beeelmanl <1986ft and iBalsara fc Krolikl (Il993ft 

suggest that the continuum from the central engine can 
cause X-ray heated winds, which contain gas that evap- 
orated off a surrounding torus. Their models show that 
the flow speeds of these winds can reach several hundred 
km s" 1 in the pc-s cale region around the AGN. Following 
iDopita et al] J2002) , also further from the nucleus the nu- 
clear radiation can be responsible, if coupled to dust, for 
radiation pressure dominated gas outflows. This idea has 
been applied to the case of NGC 1068, where it can explain 
an outflow of ~ 700 km s _1 at about 100 pc from the nu- 
cleus. A si milar mechanism has been suggested in the case of 
Cygnus A llvan Bemmel et al1l2003ft . In this case, the radia- 
tion from the powerful AGN in this radio galaxy can explain 
the outflow of ~ 170 km s _1 observed in the ionization cone 
at about 1 kpc from the nucleus. 

However, the AGN in the core of 3C 293 appears to 
be relatively weak compared with the objects considered 
above. Although obscuration certainly plays a role, we 
nevertheless find that the warm gas has a low ionization 
state in the centre of 3C 293. In particular, the [OIII]A5007 
line is weak. Alternatively, we can use the far-IR luminosity 
and assume that it is all due to re-radiated quasar-light 
instead of starlight. The far-IR luminosity of 3C 293, 
Leofim/IjQ = 10.08 iGolombek. Milev fc NeueebaueJ 



1988) , is at the lower e nd of the luminosity for quasars 
iNeugebauer et "all Il986ft . This is more than an order of 
magnitude lower than for the nucleus of Cygnus A and 
about a factor 6 lower than for the nucleus of NGC 1068. 
It seems, therefore, unlikely that the radiation from the 
nucleus of 3C 293 is strong enough to accelerate the gas 
to the large velocities observed in region E0. Also in the 
case of Cygnus A, the much faster (~1800 km s -1 ) outflow 
that has been detected along the radio axis in the [OIII] 
emission lines ( Tad hunterll99lft is not likely driven by AGN 
radiation pressure, but more likely explained as jet-induced. 

• Jet induced outflow: 

The fastest outflow of gas is seen ~1 kpc east of the 
nucleus in a region (E0) where there also is enhanced ra- 
dio contin uum due to a propa gating radio jet. At VLBI- 
scale (e.g. Be swick et al.ll2004ft this radio jet shows a dis- 
torted morphology, in cluding several bright knots, in par- 
ticular near region EO. lBeswick et al] 120041 argue that the 
radio jet in region E0 is approaching and that the contin- 
uum emission in this region might not be related to the 
fast jet itself, but rather to a low- velocity shear layer sur- 
rounding the jet, possibly created by an interaction between 
the jet and the ISM. In the central few kpc of 3C 293 
there is a large reservoir of dense, cold gas, including HI 
clouds (with individual cloud masses < few times 10 6 Mq \ 
lHaschick fc Baarj 1 19851: iBeswick. Pedlar fc Holowavl I2002T) 
as well as 1.5 x 10 10 M n of H> favans et al]ll999ft . Overall 
these results suggest that the observed fast outflow of gas is 
caused by an interaction between the propagating radio jet 
and this dense ISM. Moreover, the outflow in 3C 293 resem- 
bles cases of jet-related outfl ows in other near by powerful ra- 
dio galaxies like Cygnus A ]TadhunteJll99lll. PKS 2250-41 
JVillar-Martin et alj|l999ft. PKS 1549-79 (ITadhunter et al. I 
l200ll) , PKS 1345+12 jHolt Tadhunter fc Morganti 12003ft 
and 3C 305 iMorganti et al.ll2005ft .~ 

The total mass of warm gas involved in the outflow is 
~ 1 x 10 5 Mq, while the mass of the outflowing HI gas is 
~ 10 7 Mq. We can estimate a mass outflow-rate of the gas, 
assuming that the outflow is free-floating and spherical in 
shape around the radio knot (we assume that the outflow 
covers at least 2ir steradians as seen from the jet at the 
origin of the outflow): 

M « 2nR 2 ■ p- Fj ■ v out ■ m H , (3) 

with R the radius of the outflowing spheroid of gas, p the 
volume density and Fj the filling factor of this gas, v ou t 
the velocity of the outflowing gas and ma the mass of a 
H-atom. For the gas in region E0 that is flowing out at 
a velocity of 1000 km s _1 in a spherical shell with radius 
R = 0.75 kpc we determine M ~ 0.1 Mq yr -1 for the 
ionized gas and M ~ 23 Mq yr -1 for the HI gas. Even 
if the inner radius of the spheroid of outflowing HI is not 
larger than the ac tual size of the radio knot (ri tno t=32 pc; 
lAkuior et alJll996ft . the H I mass outflow-rate will still be of 
the order 1 Mq yr -1 . These outflow rates are conservative 
estimates, since the spin temperature of the outflowing gas 
is most likely higher than the assumed 100 K. The calcu- 
lated outflow-rates are nevertheless in agreement with the 
total mass of the outflow that we see at present. In case of 
a steady-state outflow of gas, the outflow must have been 
driven for about the past 10 6 yr in order to get a total out- 
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flow of A/hi ~ 10 7 M e and M ion ~ 10 5 M Q . During this 
time the outflowing gas has reached a distance of ~1 kpc, 
which is also in agreement with the size o f the region in 
which we detect the outflow. lAkuior et alJ <ll996t) estimate 
the age of the outer radio structure to be ~ 2.69 x 10 6 yr, 
assuming that the jet is moving freely at an advance speed of 
0.1c. Their estimate of the age of the inner radio jets (when 
due to re-started activity) is ~ 2.95 x 10 4 yr, although we 
note that there must be a large uncertainly in this estimate, 
given t he fact that the inne r radio jet is most likely rela- 
tivistic feeswick et aljl2004f) and interacting with a dense 
ISM. Nevertheless, our best estimate of 10 6 yr for the age 
of the outflow implies that the radio plasma has been in in- 
teraction with the gas for a significant fraction of the total 
lifetime of the radio source. A total mass outflow-rate of 23 
Mq yr -1 is at the low end of galaxy- s cale s tellar winds ob- 
served bv lRupke. Veilleux fe Sanders! |2002) in a sample of 
Ultra Luminous Infra-red Galaxies (ULIGs) , indicating that 
the jet-ISM interaction could be an important factor in the 
evolution of this galaxy. 



6.1.1 Details of the jet-ISM interaction 

Let us now focus on how the gas is driven out by the ra- 
dio jet, taking into account the energy-budget of the radio 
plasma. First we consider whether the radio power mea- 
sured for 3C 293 is indeed sufficient to produce the observed 
outflow. The total power of 3C 293 is Pi . 4 ghz ~ 2 x 10 25 
W Hz" 1 (Section [SJ. Following iBicknelll §002) (and ref- 
erences therein), the rate at which radiating electrons are 
supplied to the radio lobes is proportional to the energy 
flux (-Fe) of the radio plasma, so that: 



The momentu m-flux (F v ) of a jet is only a small fraction of 
its energy-flux <Bicknelll20o3) : 



3l 

Fe 



(4) 



where ki.4gh z is normally assumed or derived to be in the 
order of 10 -12 - 10" 11 for radio galaxies on the tens of kpc 
scale. Assuming Ki.4GHz = 10" 12 , and assuming that half 
of the energy flux is deposited in the eastern and the other 
half in the western jet/lobe structure, we expect that Fe is 
in the order of 10 44 erg s _1 for the eastern jet/lobe structure 
in 3C 293. From arguments mentioned in the previous Sec- 
tion, it is uncertain which part of the eastern radio jet/lobe 
structure is (and has been in the past) responsible for the 
fast outflow, and therefore also what the age of the outflow 
is. However, following the rough estimate from the previous 
Section that the outflow has been driven for about the past 
10 6 yr, the total energy supplied by the radio jet in region 
E0 over this period is in the order of a few x 10 57 erg. The 
total kinetic energy of the outflowing gas in region E0 is 
Emu = \Mv 2 w 1 x 10 56 erg. We therefore argue that the 
radio plasma in 3C 293 carries enough energy to create the 
fast outflow in region E0. 

The gas is accelerated by shocks that are created by the in- 
teraction between the jet and the surrounding medium. The 
gas outflow due to these shocks can be either momentum- 
or energy-driven. Let us first consider the momentum-driven 
case, in which the acceleration of the gas is caused by the 
impact of the working surface of the jet directly on a cloud. 
This happens when the jet is much lighter than the cloud. 
In this case only a small fraction of the jet energy flux is 
transmitted to the cloud, the rest is advected with the jet. 



F p = - x Fe (relativistic) 



F p — x Fe (non — relativistic), 



(5) 



For 3C 293 we sa w that Fe ~ 10 erg s~ for the eastern 
jet/lobe structure. iBeswick et alJ (120041) argue that the inner 
radio jet is most likely relativistic, therefore F p ~ 3 x 10 33 g 
cm s~ 2 . We can also derive the total momentum- flux of the 
outflowing gas directly from our observations: 



F p , = M xv 



6 >< 1034 (wi^r) (l 



g cm s 



(6) 



Therefore we argue that it is unlikely that the outflow is 
entirely momentum-driven. However, given the uncertainty 
in the mass-outflow rate (Section 16.11 . we cannot rule out 
this possibility. 

Alternatively, (part of) the fast outflow in region E0 could 
be energy-driven by a jet-induced lobe expansion. In this sce- 
nario the gas is swept up and compres sed as the radio jets 
hollow out a cocoon-lik e structure (e.g. 
iTadhunter et al. 1 1200 ill . According to 
Jl999h outflow velocities > 1000 km s can more easily be 
explained if the clouds are entrained in hot, shocked gas that 
expands out behind the bow-shock, simi lar to the expan- 
sion of an interstellar bubble (see also e.g.JStOTie^^^OTman| 



Capetti et aHll999l: 



Villar-Martin et alJ 



1992llDai fe Woodwardll994l:lKlein. McKee fe Colellalll994l7 



O'Dea et alj|2003h . This process has alr eady been suggested 
for several other radio galaxies (e.g. IVillar-Martin et alJ 
l!999HBest. Rottgering fe Longairl200q) . Such a more gentle 
energy-driven mechanism might also be better for explaining 
the large amounts of neutral gas involved in the outflow. 



6.1.2 Neutral gas involved in the outflow 

The fact that only about 1 per cent of the outflowing gas 
appears to be ionized, while the rest is in its neutral state, 
is maybe one of the most intriguing results from this paper 
and Paper 1. The question is how, despite the high energies 
that must be involved in a jet-ISM interaction, such large 
amounts of gas stay, or become again, neutral? 

The answer may lie within various simulations of 
jet-cloud interactions. For exam ple, simulations by 
(Sutherland. Bicknell fe Dopital i2003T) show that post-shock 
gas overtaken by a radio jet can show complex cooling. 
In the fractal structure of the gas dense fila ments can be 
forme d. Takin g into account cooling effects , ISaxton et alJ 
J2005I) and iBicknell. Saxton fe Sutherland! (120031) show 
that in various jet-ISM interactions intervening clouds can 
severely disrupt the radio jet, while at the same time some 
clouds are accelerated away from their initial positions 
by the thrust of the radio jet. Large clouds, in particular 
the ones not directly in the path of the jet, can survive 
for a long time after the rad i o jet p assed. Simulations by 
iMellema. Kurk fe Rottgerind i2002T) investigate the fate 
of a cloud of ionized gas that is overrun by the lobe of 
an expanding radio jet. Due to an underpressure in the 
cloud compared to the overpressured cocoon, shocks start 
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to travel into the cloud and the cloud is compressed and 
fragmented. Due to rapid cooling the bulk of gas gets 
locked up in these dense fragments that can survive for a 
long time and reach speeds up to 500 km s -1 . The further 
evolution of these fragments will be dominated by further 
acceleration and erosion by the passing flow, as well as 
gravitational collap se (which ul timately will result in star 
formation). iFraeile et al.l i2004) reach similar conclusions 
from their models on jet-induced star formation. They show 
that in the interaction of shocks with a radiative cloud, a 
large fraction of the cloud-gas cools to low temperatures. 
The neutral hydrogen gas in the starforming Minkowski's 
Object, for example, may have cooled from the initial warm 
gas as a result of the radiative cooling triggered by a passing 
radio jet. The amount of outflowing neutral hydrogen gas 
that we find in region E0 in 3C 293 might be explained in 
a simi lar way, although in the simulations bv lFraeile et al] 
( 2004) the post-shock gas is only accelerated very slowly. 

In light of the fast outflow of neutral hydrogen gas that we 
find in 3C 293, it would be interesting to explore whether, 
in the above mentioned simulations, at least a significant 
part of the dense clouds or fragments can consist of gas 
that remained or recombined again into neutral hydrogen 
gas, and if this gas can reach velocities up to 1000 km s" 1 
during the lifetime of the radio source. 



the much fainter large-scale jet (iMorganti et al.lll999h . In 
the case of Centaurus A it has been proposed that the large 
velocities shown by the gas are still created by the strong 
instabilities produced by the propagation of the jet in the 
past (especially in the transition region between the inner 
and the outer lobe). A similar scenario could apply to the 
outflow observed in region E2 in 3C 293. 

Using similar arguments to those used for region E0, 
the outflow in region E2, although it appears to have bro- 
ken through the ionized disk of the galaxy, is still not 
spread widely enough that it can be explained by a fossil 
wind from the post-st arburst stellar population found by 
Tadhuntcr ct al Als o, from both our op tical spectra 

as well as the UV-data from lAllen et all i2002l) . there is no 
indication that a substantially younger stellar population or 
HII region is present in region E2. Therefore it is not likely 
that the outflow in region E2 is driven by stellar winds. 

An AGN-wind is not likely strong enough to explain 
the outflow in region E2, since the outflow seems to be too 
extreme for that. However, for the outflow in region El and 
in particular regions Wl and W2 (assuming the outflow is 
real in these western regions and not a result of observational 
line-of-sight effects as discussed in Section 13. H we cannot 
discard an AGN-wind as the possible driving mechanism 
with the available data. 



6.2 Outflows in other regions 

As explained in Section ^. II there are other regions in which 
complex kinematics of the gas, possibly connected to out- 
flows, are seen. However, the gas kinematics observed in the 
other regions are not as extreme as in region E0. 

The propagating radio jet, responsible for the fast out- 
flow in region E0 as explained above, likely carries enough 
energy to also cause outflows in the regions other than E0. 
This is an appealing scenario, in particular since the com- 
plex kinematics are observed along the inner radio axis. 
The fact that the broad component of the emission-line gas 
(the component that represents the ionized gas outflows) is 
blueshifted w.r.t. the narrow component on the eastern side 
and redshifted on the western side (in case the asymmetric 
emission-line profile in regions Wl and W2 represents a real 
outflow and is not the effect of observational line-of-sight ef- 
fects as discussed in Section ^. 11 , implies that the radio jet is 
approaching on the eastern side and receiding on the western 
side. However, in region E2 so far no radio continuum has 
been detected (see Section f3.lL although it is the region in 
between the bright inner radio jet and the bright outer radio 
lobe. This could mean that the radio jet is propagating more 
efficiently in this region and it might suffer from Doppler 
de-boosting. Something similar likely happens in the region 
between the core and the region E0 as observed with very 
high resolution VLBI imaging feeswick et al.ll2004l) . An- 
other possibility is that the inner jet repre sent re-started 
activity (as suggested bv lAkuior et al.l [l996'). and the radio 
continuum in between the inner jet and the outer lobe has 
faded away. This last explanation could resemble the case 
of Centaurus A. As in 3C 293, also in Centaurus A a so 
called i nner-filament of ionize d gas with very complex kine- 
matics l|Morganti et al.l fl991'l is located outside the inner 
bright radio lobe (the lobe that can be considered to be the 
result of the most recent activity), but about 2 kpc from 
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